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We study the spectrum of fermion states localized within the vortex core of a weak-coupling p- 
wave superfluid. The low energy spectrum consists of two anomalous branches that generate large 
density of states at the locations of the half cores of the vortex. Fermi liquid interactions significantly 
stretch the vortex structure, which leads to Lifshitz transition in the effective Fermi surface of the 
vortex core fermions. We apply the results to rotational dynamics of vortices in superfluid 3 He-B 
and find explanation for the observed slow mode. 


PACS numbers: 67.30.he, 67.30.hj, 74.25.nj 

The double-core vortex is an amazing structure be¬ 
cause it is the unique answer to a simple question: what 
is the vortex structure of a weak-coupling p-wave-pairing 
superfluid. The ground state in this case is the Balian- 
Werthamer (BW) state [lj • The lowest energy vortex has 
the double-core structure, where the core is split into two 
“half cores” as depicted in Fig. HK a ,b) M. This is not 
only of theoretical interest since superfluid 3 He is close 
to being weak coupling, and its B phase was identified 
as the BW state. Two vortex types have been found ex¬ 
perimentally in 3 He-B [Ml- The vortex being stable in 
the major, low-pressure part of the phase diagram was 
identified as the double-core vortex. The vortex stable at 
higher pressures was identified as the axially symmetric 
A-phase-core vortex. Available experimental evidence is 
consistent with the theoretical identification of the vor¬ 
tex structures. In particular, the broken axial symmetry 
of the double-core vortex was used to explain the pe¬ 
culiar dynamical properties that have been observed for 
the low pressure vortex using homogeneous precessing 
domain (HPD) mode of NMR ul flij . A similar double¬ 
core vortex structure has been suggested to appear in 
spin-triplet heavy fermion superconductor UPt 3 |15| . 

One of the most interesting properties of quantized 
vortices in superconductors and Fermi superfluids is the 
presence of fermionic quasiparticles localized within vor¬ 
tex core at energies smaller than the bulk energy gap 
0, S3- Generally fermionic bound states determine 
both thermodynamic and dynamic properties of vortices 
at low temperatures jl8l-l22j. In the rotational dynam¬ 
ics of the double-core vortex they are predicted to give 
rise to resonance absorption at the frequency comparable 
with the spacing of the localized energy eigenstates [23| . 
Recently much attention has been focused on the topo¬ 
logically protected zero energy vortex-core and surface 
states in superfluid 3 He [24M27| . Particularly motivat¬ 
ing is a predicted existence of self-conjugated Majorana 
states localized on half-quantum vortices in p-wave su¬ 
perfluids [H]. 

In this letter we calculate the low-energy fermionic ex¬ 
citation spectrum of the double-core vortex. We find that 
the low-energy excitations mostly are localized in the two 
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FIG. 1: (Color online) (a,b) The double-core vortex structure 
made visible by the pair density I'Ll 2 = i |A M ,i| 2 plotted in 
the x—y plane perpendicular to the vortex axis at temperature 
T/T c = 0.9, and (b) T/T c = 0.1. (c,d) The normalized local 
density of states profiles demonstrating the quasiparticle wave 
function at the Fermi level e = 0, p z = 0 at T/T c = 0.9 and 
(d) T/T c — 0.1. All plots correspond to pressure P = 24 bar. 


half cores. This is visualized in Fig.[l] c,d) which show the 
fermionic local density of states (LDOS) profiles around 
the vortex core. We can interpret the two half cores 
as potential wells for quasiparticles. The motion of the 
excitations between the wells depends on the potential 
barrier between them. We find that this barrier changes 
essentially as the distance of the wells changes as a func¬ 
tion of pressure and temperature (see Fig. |T] to compare 
vortex structures at T/T c = 0.9 and T/T c = 0.1). This 
implies a transition from excitations that circle both half 
cores to separate excitations that circle only a single half 
core. This transition can be seen as a Lifshitz-type tran¬ 
sition in the topology of the fermionic states bound to 
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the vortex core. We discuss how this could be observable 
in the time scales of rotational dynamics. Comparing our 
calculations with an earlier experiment reveals a serious 
disagreement in the model used to interpret the experi¬ 
mental data 0. We construct a different model, which 
also provides an explanation for the l ong time scale ob¬ 
served in rotational dynamics 0,0, [30]. 

The triplet pairing of fermions in orbital p -wave states 
is described by the matrix 


A (r,p) = A ai (r)i& a a y pi , (1) 

OL,i 

where b X)y<z are Pauli matrices, p is the momentum close 
to the Fermi surface p ps pf = hkp, and p = p/p. The 
gap function A CQ is determined by the 3x3 order 
parameter matrix with complex components A a ;. Here 
a = x,y,z and i = x,y, z are spin and orbital indices, 
respectively. 

In the weak coupling theory of p- wave superfluid, the 
stable state has the Balian-Werthamer (BW) form JJ. 
In BW state, the order parameter far from the vortex 
axis is A a i = Ao exp (icp)R a i. Here tp is the azimuth with 
respect to the vortex axis, R a i is a constant rotation 
matrix and Ao the order parameter amplitude. Near the 
vortex axis a more sophisticated structure appears @- 
[s[. It is energetically favorable to change the sign of the 
order parameter across the vortex axis by spin rotation 
of the BW-state matrix A a i by tt Q. This effectively 
results in splitting of a singly-quantized vortex to a pair 
of half-quantum vortices that are bound together by a 
planar-phase domain wall. For illustration see Figs.[L|a) 
and[l]]b), which show the pair density |T| 2 = J2 a i |4lai| 2 
in the x — y plane. The pair density has two distinct 
minima, whence the name double-core vortex. 

To determine the vortex structure we calculate self- 
consistently the order parameter and the Fermi-liquid self 
energy [311 ]. The numerics is performed as described in 
Ref. [a, i.e. using the explosion trick to solve the Eilen- 
berger transport equation. We extend the previous work 
@,:7] to higher accuracy, lower temperatures and different 
values of the Fermi-liquid parameter F{ corresponding 
to different pressures. The parameter Ff determines the 
feedback of superfluid mass current on the order parame¬ 
ter and can significantly change both the vortex structure 
and spectrum of bound fermions. 

The distance a between the half cores is shown Fig. [2] 
Its scale is Rq = (1 + .Ff/3) £o, where £o = h,VF/^T c is 
the coherence length and vf the Fermi velocity. As F/ 
in liquid 3 He ranges from 5.4 to 14.6 depending on pres¬ 
sure P [32|], the two length scales can differ essentially. 
Thus at large values of corresponding to high pres¬ 
sures, the vortex size at a low temperature is much larger 
than the coherence length. For example, a = 46£o in the 
case of Fig. QJb). Fig. [3] also shows strong temperature 
and pressure dependence. The distance of the half cores 
grows almost 3-fold when the temperature decreases from 
0.9T C to 0.1T C at 24 bar. Similarly as a function of pres¬ 
sure the distance a measured in units of £ 0 , grows almost 
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FIG. 2: (Color online) Distance a between the half cores in 
the double-core vortex at different values of the Fermi liquid 
parameter F*. The locations of the half cores are determined 
from zeros of supercurrent density. 


2-fold when the pressure increases from 0 to 24 bar at a 
low temperature. 

We calculate the excitation spectrum of the double¬ 
core vortex using the self-consistent order parameter 
field. This is done by solving the eigenvalue problem 
for the system of Andreev equations, which are ordinary 
differential equations describing the propagation of quasi¬ 
particle wave function along classical trajectories. 

The momentum p of a low energy excitation is close 
to the Fermi surface, p ~ pf- The classical trajecto¬ 
ries are straight lines parallel to p. In studying a vortex 
we fix the z axis as the vortex axis, and we parameter¬ 
ize the momentum direction p = (pj_ cos 0 P ,p± sin 0 p ,p z ). 
The direction on the trajectory is fixed by giving p z and 
Op. The location of the trajectory is given by the im¬ 
pact parameter b , the coordinate measuring the distance 
from the vortex axis. The parameterization is visual¬ 
ized in Fig. [U]a). The impact parameter is related to 
the projection of the angular momentum p on the vor¬ 
tex axis through the usual classical mechanics formula 
p = p±b. The quasiclassical energy spectrum is given 
by e = £i(p z , 0 p ,b), where the parameters p z ,0 p ,b specify 
the classical trajectory and integer i counts the eigen¬ 
values of Andreev equations on a given trajectory [331 ]. 
Figure (Hb) shows a bunch of trajectories at the Fermi 
level and p z = 0. The concentration of the trajectories 
at the two half cores results in the large LDOS at the 
half cores. The concave triangular shape of the caustic 
of the trajectories at the half cores is clearly visible in 
the LDOS shown in Fig. [Hd). Also the classically non- 
allowed region around the vortex axis in Fig. [3jb) can be 
recognized in Fig. [Ijc) as a valley in the LDOS profile in 
the region between the half cores. 

Due to the lifted spin degeneracy, singly-quantized vor¬ 
tices in 3 He-B have two anomalous branches of quasipar¬ 
ticle spectrum Q. At low energy compared to the bulk 
energy gap, |e| <C Ao, they can be represented as 

£i ( p z , Op, b) = -uJiPF ( b-bi ), (2) 

where i = 1,2. Here bi(j> Z: 0 p ) is the impact parame- 
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FIG. 3: (Color online) (a) Schematic plot of a quasiclassical 
trajectory in the x — y plane in the direction of (cos 9 P , sin 9 P ) 
passing the vortex axis at distance b (impact parameter). A 
point r = (x , y ) on the trajectory is determined by the coor¬ 
dinate s. (b) A bunch of quasiparticle trajectories (straight 
lines) at the Fermi level e = 0 and p z — 0 superimposed 
on the pair-density contour plot at T/T c = 0.5 and F‘ = 0. 
(c) Anomalous branches of the quasiparticle spectrum e = 
ei ,2 (b,9 p ) at p z = 0, F" = 12 and T = 0.9T C . (d) Two sheets 
of effective Fermi surface b = 61,2 (9 P ) at p z = 0, F( = 12, 
T = 0.9T C (solid lines) and T = 0.1T C (dashed lines). 


ter that corresponds to vanishing excitation energy, and 
Ui(p z ,0p) indicates the slope of the energy at b = b,. The 
new feature in a non-axisymmetric vortex is that these 
parameters depend on the trajectory direction 9 p in the 
x — y plane. 

Fig- He) shows the calculated quasiparticle energies as 
a function of impact parameter b and different directions 
of the trajectory 9 p . The curves cross the Fermi level at a 
finite b in accordance with Eq. ©• These locations in the 
energy spectrum bi^{9 p ) are shown by the black dots in 
Fig. He). The states at the Fermi level in the spectrum 
© form a 2D effective Fermi surface b = 61 , 2 {Pzi & P ) 
in the 3D space formed by the quasiclassical quantum 
numbers ( p z ,9 p ,b) in the vortex core. Because of two 
nondegenerate branches ©, there are two sheets in the 
Fermi surface. One more representation of this is given 
in Fig.Hd)- It shows as a function of 9 p . The curves 
depend also on p z but that dependence is less important 
in the following because p z is conserved. For comparison, 
the trajectories passing precisely through the half cores at 
y = ±a/2 would correspond to curves b(9 p ) = =F \a cos 9 p . 


The topology of the effective Fermi surface is deter¬ 
mined by the behavior of zero energy lines £q, 2 ($p) at 
9 P = tt( n+ 1/2) with n = 1,2. At these angles the quasi¬ 
particle trajectories pass through both half cores. In gen¬ 
eral there is overlap of the quasiparticle wave functions 
localized at different half cores. This makes that there is 
no sign change of 61,2 (^fo)- That is, there is anticrossing of 
the two branches and a finite splitting 25b = \ b\ — 62 ) > 0 
at 9 p = n/2, as shown by the solid lines in Fig. Hd). 
Physically this means that an excitation created at one 
half core will jump periodically between the half cores. 

The growing core separation (compared to £ 0 ) a t I° w 
temperatures and large pressures reduces the overlap of 
the quasiparticle wave functions located at different half 
cores. As a result the splitting 25b becomes extremely 
small as shown by dashed lines in Fig. Hd) for Ff = 12 
and T = 0.1T C . In this case Landau-Zener (LZ) tunnel¬ 
ing between the quasiclassical branches © becomes im¬ 
portant. The probability W of these transitions can be 
found from the conventional approach [3|J, [36| by taking 
9 p and the angular momentum p = p±b as the conju¬ 
gate variables. Near the anticrossing point at 9 p = 7t/2 
we can approximate 61,2 (0 P ) « ±\j5b 2 + (a9/ 2) 2 , where 
9 = 9 p — 7t/2. The transition probability is given by 

W = exp[—2fcj_Im f^ 8 ( b± — 62 )^] where i9* = 2 i5b/a 
is the intersection point of quasiclassical branches in 
the complex plane. A simple calc ulation yields W = 
exp[— 2irp±(5b/ Ab) 2 }, where A b = ffa/kF has a physical 
meaning of the quantum mechanical uncertainty of the 
impact parameter. 

Once the transition probability becomes large, W ss 1, 
the LZ tunneling changes the topology of the effective 
Fermi surface so that quasiparticles remain localized in 
one or the other of the half cores. In Fig. HA) this means 
transition to intersecting zero energy curves £> 1 ( 2 ) ( 0 p) = 

£> 1(2 )(0 P ) for -7t/2 < 9 p < it 12 and 61 ( 2 ) (0p) = £> 2 ( 1 ) (#p) 
for 7t/2 < 9 p < 3-7t/2. The calculated LZ probability 
W(T, P,p z ) is shown in supplemental material [37j to 
demonstrate that the condition W ~ 1 is realized in the 
double-core vortex at large pressures and low temper¬ 
atures. The crossover from split b\^{9 p ) to intersecting 
isoenergetic lines £> 1,2 (9 P ) is an analog of the Lifshitz tran¬ 
sition |39| changing the topology of the Fermi surface. 

The transition leads to a formation of two spatially sep¬ 
arated low-energy fermionic states localized at the half¬ 
cores. Whether there are Majorana states precisely at the 
Fermi level |28| or not H3] is beyond our quasiclassical 
approach. The transition affects the rotational dynamics 
of the double-core vortex. The bound fermions in the core 
respond to oscillation of the core orientation. A friction 
torque acting on a rotating vortex core can be expressed 
by a friction coefficient / = /iPf(^fCo) 2 , where /1 ~ 1 
is dimensionless and the factor pf(^fCo ) 2 is determined 
by the density of quasiparticles in the vortex core. The 
expression for the friction torque (23| yields resonance 
peaks in / located at angular frequencies ui « nE m /H 
where n is integer. Here E m is the spacing of quantized 
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P = 2 bar 

P = 12 bar 

P = 24 bar 

T = 0.05T C 

27 kHz 

71 kHz 

98 kHz 

T = 0.5 T c 

22 kHz 

65 kHz 

106 kHz 


TABLE I: Values of the minigap Em/h at different pressures 
and temperatures. 




FIG. 4: (Color online) Demonstration of the effect of the Lif- 
shitz transition on the rotational friction coefficient /i plotted 
as a function of frequency oj for (a) W = 0, (b) W = 1. The 
vortex structure and the excitation spectrum are calculated 
at T = 0.5 T c and (a) P = 2 bar, (b) P = 24 bar. The mini¬ 
gap values are given in Tabic [T] According to mutual friction 
measurements [12( h/E m T = 0.7 at P = 24 bar but a wider 
range is given to illustrate the influence of relaxation on the 
shape of absorption peak. 


energy levels obtained from the quasiclassical spectrum 
Q using the Bohr-Sommerfeld quantization rule for the 
angular momentum [23[, E m = h{uj^ 1 {p z = O)}^ 1 , where 
(...) denotes the average over 9 p . The scale of the minigap 
is determined by h/T n = (2ttT c ) 2 /vfPf, which is on the 
order of the quasiparticle relaxation rate in the normal 
state. The calculated values of the minigap are listed in 
Table U 

The amplitudes of the resonances are determined by 
the Fourier amplitudes of the zero-energy curves shown 
in Fig. [3] A n ~ | f dO p e m6p bi(9 p )\. From the plots in 
Fig. 0(d) one can see that at pressures below the Lifshitz 
transition the largest components are those with double 
frequency Huj = 2 E m . At pressures above the transition 
the amplitudes are determined by the harmonics of the 
intersecting curves bi{9 p ) : which have strongest matrix 
element at Hlo = E m . The difference in the friction co¬ 
efficient fi in the two cases is demonstrated in Fig. [4] 
In the following we show that at least the low frequency 
limit of the curves in Fig. 0] is experimentally accessible. 

Kondo et al [H, HH have studied a sample of rotating 
3 He-B using the homogeneously precessing domain. In 
this mode the magnetization M is tipped by a large an¬ 
gle (> 104°) from the held direction B. It was found that 
the contribution of vortices to the relaxation changed on 
a few minute time scale pFjj]. This was interpreted that 


the double-core vortex gets twisted as its end points (at 
2 = EL/2) are pinned but in the bulk the rotating mag¬ 
netization exerts a torque on the core. A quantitative 
model was constructed for the vortex core rotation angle 
z ) as a function of time t and 2 . The parameters of 
the model were determined by fitting to the experiment 
fl4| . These include the friction parameter /, the dipole 
torque Tf >, which drives the vortex in the presence of ro¬ 
tating magnetization, and the rigidity AT, which gives the 
energy caused by twisting the core, Ftwist = \K{d z (j)) 2 . 

By precise calculation of the vortex structure we can 
now calculate the vortex parameters. We find a value of 
/ that is three orders of magnitude larger than fitted by 
Kondo et al 0 - Thus a serious revision of the model 
has to be made. The large value of / means that only a 
negligible fraction of energy dissipation comes from vor¬ 
tex core rotation. Thus essentially all dissipation has to 
arise from normal-superfluid disequilibrium [40l | , spin dif¬ 
fusion, and radiation of spin waves. Without going into 
details, these can be incorporated by allowing an elastic 
vortex structure, where the rotation angle a(t, z) at a dis¬ 
tance from the vortex axis (where the dipole torque acts) 
can be different from the vortex core angle </>(t, z). These 
are bound by elasticity energy ^Ta(ol — </>) 2 , and both an¬ 
gles have their own friction coefficients: f<p = — SF/Scj. i, 
ga = —6F/6a. This model results in the diffusion equa¬ 
tion Iii7| 


jd 2 z cfr+^-, 

f vf 


(3) 


where P g is the power absorption per vortex length. An 
important virtue of this model is that based on our cal¬ 
culations of / and K, Eq. 0 predicts the time scale 
L 2 f/n 2 K of several minutes. Thus Eq. ([3]) give s a sim¬ 
ple explanation for the observed slow mode [ll| 29, 
which remained unexplained in previous models [14, 


In summary, we have investigated the spectrum of 
bound fermion states localized within the vortex core 
of weak-coupling p-wave superfluid. We predicted a Lif¬ 
shitz transition, which separates low-energy quasiparticle 
states at the half cores and affects the rotational dynam¬ 
ics. Applying our results to Ref. [3 explains the observed 
long time scale and thus gives one more piece of evidence 
of the double-core nature of the low pressure vortex in 
3 He-B. 
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I. LANDAU-ZENER TUNNELING PROBABILITY BETWEEN THE QUASICLASSICAL SPECTRUM 

BRANCHES 

The topology of quasiparticle spectrum in double-core vortices is determined by the behavior of quasiclassical 
spectrum branches e = eqRA, 9 p ,b) near the intersection points at 9 P = 7 t /2 + 7 rn, see Fig. 3(d) in the main text. 
In the Fig. [5] we show in detail these curves in the vicinity of an anticrossing point 9 p = 7 r /2 for several values 
of parameters. This plot illustrates main properties of the splitting 5b as function of temperature T, pressure P, 
momentum projection to the vortex axis p z . 

Quantitatively the splitting 5b is determined by the overlap between quasiparticle wave functions localized at 
different half-cores as shown in the Fig. 1 in the main text. At low temperatures the characteristic localization scale 
is determined by the coherence length £o- On the other hand the distance between vortex cores a is determined 
by the scale Po = (1 + Pf)£o so that the ratio I?o/£o increases with growing pressure. Hence the overlap at large 
pressures is weaker and the splitting 5b decreases as can be seen comparing the curves for Pj 3 = 9 (P = 11.6 bar) 
and F* = 12 (P = 23.75 bar) in Fig. ER)- At larger temperatures the quasiparticlc localization is determined by 
the temperature-dependent coherence length £ s=s Hvf/ A where A = A (T) is the gap amplitude. On the other hand 
the distance between vortex cores for temperatures up to T = 0.99T C has a much weaker temperature dependence for 
Pf = 6 , 9, 12 as shown in Fig. 2 in the main text. Thus the overlap between half-core states and the splitting 5b 
strongly decrease with decreasing temperature which can be seen from the comparison of the curves b = b\ 2 (9 V ) for 
T = 0.9 T c and 0.1T C in Fig.EIM . 

The most striking is the dependence of splitting 5b on the quasiparticle momentum projection on the vortex axis 
p z . The absolute value of momentum is fixed and determined by the Fermi momentum. Hence different values of p z 
correspond to the different angles of quasiclassical trajectories with respect to the vortex axis. For the finite p z the 
effective distance between half-cores along the trajectories is elongated by the factor of l/-y/l - p%. which decreases 
the overlap of localized states and suppresses the splitting 5b. As can be seen in Figs.EU) and Eld), the splitting can 
decrease to the very small values at large p z ~ 1 . 

The splitting 5b determines the tunneling between the quasiclassical branches. The Landau-Zener tunneling prob¬ 
ability W = W(T,P,p z ) calculated according to the general expression from the main text is shown in Figs. El a) and 
Elb) as function of temperature for several values of P and p z . These curves demonstrate a general tendency of the 
tunneling probability to increase with decreasing temperatures and increasing pressure. The probability is strongly 
enhanced at large p z . Thus the Lifshitz transition determined by W ~ 1 occurs at different values of P and T for 
different p z . From Fig. El one can see that at high pressures and low temperatures the condition W ~ 1 is valid for all 
values of p z . 


II. MODEL OF ROTATIONAL DYNAMICS 


The rotational dynamics of the double-core vortex was studied by Kondo et al HI. Here we present a modified 
model that allows azimuthal shear of the vortex structure. Such a case appears because the driving dipole torque 

acts on the asymptotic order parameter, typically a distance £d ~ 10 pro. from the vortex axis, whereas the opposing 

rotational friction occurs in the vortex cores, on the scale of coherence length £0 ~ 10 nm. Our notation follows closely 
Ref. 0 and the model of Ref. □ is obtained in the limiting case Ta — > 00 and g = 0 . 

We study the model where the free energy of a single vortex parallel to z is 


F = 


dz 


—T H {S a R a ibi) 2 + T D (a ■ hj 2 


T A {d ■ b ) 2 + ^K{d z jg) 2 


(4) 


The magnetic field term (coefficient Th ) depends on the hard core anisotropy axis b = x cos <f>+y sin </>. The dipole term 
(coefficient To) depends on the soft core anisotropy axis a = a; cos a + y sin a. The azimuthal shear term (coefficient 
T a ) is supposed to be strong enough to keep the two anisotropy vectors nearly parallel. Other quantities are the spin- 
orbit rotation matrix R a i parametrized by axis h = y cos cot + (z sin y — xcosrj) sin tut and angle 9 = arccos(—1/4), 
the direction of the precessing magnetization S a = R a iHi , and the static field H = F[(xsinr] + zcosy). This model 
neglects the anisotropy of the dipole energy in the double-core vortex in the plane perpendicular to a s=s b. In principle, 
we should allow a to have component in the z direction also, but this will not affect the main results and thus is 
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FIG. 5: (Color online) Detailed structure of zero-energy lines b = bi, 2 (d p ) near the anticrossing point 6 P = n/2 for different 
values of parameters, (a) T = 0.1 T c and different values of pressure Ff — 9 (red), FI = 12 (blue); (b) FI = 12 and T = 0.9 T c 
(blue), T = 0.1 T c (red); (c) FI = 12, (d) FI = 9 and different momentum projections p z = 0 (red), 0.23 (blue), 0.45 (green). 


dropped here for simplicity. We get 


= dz 


F = J dz —T h {H ■ b ) 2 + T D (a ■ n ) 2 - T a (a ■ bf + - K(d z 0) 2 
—Th sin 2 pcos 2 cj) + Td (cos ut sin a — cos 77 sin wf cos a) 2 + T A {a — </>) 2 + ^-K(d z (j)) 2 


(5) 


Note that for H = x, cos 77 = 0 and both field and dipole terms are minimized by <f> = a = nit, that is a = b = ±x 
and h rotating around it. 

We suppose the frictional equations of motion 


SF ■ 6F 

oa oq) 


( 6 ) 


We obtain 


SF 


ga = — ^r— = —Td sin 2 uit sin 2 77 sin 2a — Td cos2wf(l — cos 77 ) sin 2a + Td cos 77 sin 2 (cot — a) — 2 T A (a — </)). (7) 

oa 

Using the same approximation as Kondo et al this simplifies to 

ga = s ^ n2 r ? s i n 2 o: + Td cos 77 sin 2 (wt — a) — 2 T A {a — (f>). 


( 8 ) 
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FIG. 6: (Color online) Temperature dependencies of the Landau-Zener tunneling probability between quasiclassical spectral 
branches in the double-core vortex for (a) Ff = 9 (P = 11.6 bar), (b) FI = 12 {P = 23.75 bar) and different values of 
momentum projection p z . 


This expression is exact in second (but not fourth) power in 77 . The other equation is 

8F 

f(j) = —— = —T h sin 2 rj sin 20 + 2 T A (a — 0) + A"<9 2 0. (9) 

dcp 

Suppose now that / is very large, which means (f>(t) is so slow that it can be taken as a constant during one cycle. 
We study 

ga = —M sin 2a + Gsm2(ojt — a) — 2 T^{a — <j>) (10) 

where in the approximation above M = \Td sin 2 g and G = Td cos 77 . Suppose a(t) = «o + ou(t) with small ctq 

ga 1 = — M sin2a 0 ~ 2Tyi(ao — 4>) + Gsin2(u;t — ao) — 2 (Ta + M cos2ao)ai — 2Gcos2(wt — olo)ol\. (11) 
Assuming a\ = A cos 2 (cut — \) 

— 2ujgAsin2(ujt — \) 

= —M sin2a 0 _ 2Tyi(ao — 4>) + Gsin2(wi — ao) — 2 A{Ta + M cos2ao) cos2(wt — x) 

—GAcos2(2wt — ao — x) ~ GAcos2(y — ao). (12) 

At frequencies 0 and 2 u> we get 


0 = —M sin2ao — 2TA{ao — </>)— GAcos2(y — ao) 

—2ujgAsm2(ujt — x) = Gsin2(wt — ao) — 2 A(Ta + M cos2ao) cos 2 (cot — y) (13) 

The latter equation can be written 

2A(Ta + M cos 2ao) cos 2 (cot — y) — 2ujgA sin 2(cot — y) = G sin 2{ujt — ao) (14) 

2A(Ta + M cos2ao) cos 2 wt — 2ujgAsm2ujt = Gsin 2 (wt + y — ao) 

= G[sin2wicos2(y — ao) + cos2wtsin2(y — ao)] (15) 

2 A{Ta + M cos 2ao) = G sin 2(y — ao) 

— 2ujgA = Gcos2(y — ao) (16) 
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A 2 = 


G 2 /4 


G 2 /4 _ lT^cos 2 r] 


u> 2 g 2 + ( Ta + M cos 2ao) 2 to 2 g 2 + T\ uj 2 g 2 + T\ 

The absorbed power is determined by A, 

p g = = 9(a 2 ) = 2 gu 2 A 2 . 

Note that here P denotes the absorption per vortex length whereas in Ref. [3 P denotes the total absorption. 
Applying the latter of (flGl) to the first equation (fl3l) gives 


From (O we have 


M sin 2ag + 2Ta(o-o — 0) — 2ujgA 2 = 0. 


/0 = —T h sin 2 77 sin 20 + 2Ta(u 0 — 0 ) + Kd 2 (p. 


(17) 


(18) 


(19) 


( 20 ) 


Because 0 is a slow variable, it is only the low frequency limit of / that appears. What remains to be solved is 0 and 
ao from equations uni) and Two different cases need to be studied. 

1) Rocking oscillations. We assume a solution where the left hand side of (l20l) vanishes, 


SF 

—jTT- = —T h sin 2 77 sin 20 + 2Ta(cxo — 0 ) + Kd 2 (j) = 0 . 
00 


Supposing 0 - a <1 we can solve this 


1 = ao 


Th sin 2 ? 7 sin 2 ao — Kd 2 a 0 


ao 


Th sin 2 77 sin 2 a 0 — Kd 2 a 0 


2 (Th sin 2 77 cos 2ao + Ta) u 2Ta 

The lag of 0 behind oq is an increasing function of 0 ps oq for 0 < 0 « ao < 7t/ 4. Substituting 


(Th + -To) sin 2 77 sin 2a 0 — Kd 2 oto — 2u>gA 2 = 0 


and thus 


sin 2 ao = 


2 cog A 2 + Kd 2 ao 
(T h + \T d ) sin 2 77 


( 21 ) 

( 22 ) 

into (flop gives 

(23) 

(24) 


This solution is possible only if the right hand side of (l24l) is less than unity. At equality, ao = 7 r /4 + 7 i 7 r. The stability 
condition at 9 2 ao = 0 can be written 


tan 2 77 > tan 2 770 = 


\u9Tl 


(T H + \T D )(u 2 g 2 + T 2 A y 


Interestingly, this can also be written 


sin 2 77 > sin 2 ? 7 0 = 


Pn 


u(T h + \t d ) 


(25) 


(26) 


These results are essentially (tan 770 ~ sin 770 ) the same as in Ref. [3 if on e makes the replacement / —> g + T A /gu> 2 . 

2) In the case solution (l24l) is not possible, we get slow rolling motion of 0(£). For simplicity we consider 77 = 0 
only. From (fl9)l we get 


and substitution to m gives 


Equation 0 is a diffusion equation 


2 T A (a 0 - 0 ) = 2 ujgA 2 


-7 *+& 


(27) 


© 


0 = Dd 2 z (j) + c 


( 28 ) 
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The solution for 0 < z < L and t > 0 with boundary conditions <j>(z = 0) = 0 = <p(z = L) can be found as Fourier 
series 


& = z 2 ) + ^2,B n e “"‘sin k n z (29) 

n 

with k n = n7r/ L. The slowest component has the rate 

T = tt 2 D/L 2 = n 2 K/fL 2 . (30) 

We compare our results to Ref. [Til . For that we give numerical values of our weak-coupling results corresponding 
to 29.3 bars pressure, T/T c = 0.5 and L = 7 mm. At this pressure Ff = 13.3 and we neglect other Fermi-liquid 
interactions. For the friction coefficient we get / w 2 x 10 -19 J s/m, see Section Hill In section HVl we calculate 
K « 5 x 10” 27 Jm. This gives the time constant T _1 s=s 3 min. This order of magnitude is in good agreement with 
the measured value of the “slow mode” in Ref. [l3|, [29|, [33- Extensive measurements of the slow mode are presented in 
Ref. M It seems that these results need to be reinterpreted using the diffusion equation ([3]). 

Our results could be compared to / = 4.5 x 10 -23 J s/m and K = 1.1 x 10~ 25 Jm that were obtained by fitting 
the experiment to the model in Ref. [lj- 


III. ROTATIONAL FRICTION 


The calculation of the rotational friction coefficient / in the main text applies the kinetic equation of Ref. |23| to 
our numerical solution of the excitation spectrum in the vortex core. The results are presented in Fig. [4] Evaluating 
this at the conditions of Ref. Hi (see above) gives / = 1.7 x 10 19 J s/m in the low frequency limit. The result at the 
Larmor precession frequency at H = 14.2 mT ( huj/E m = 3.9) would be /(wi) = 0.9 x 10“ 19 J s/m. 

An alternative estimate of / is obtained as follows. In the case of strong Landau-Zener tunneling, the half cores 
behave as separated half-quantum vortices except that they are bound together at the distance a. The standard 
mutual friction force [HI, [HI applied to both half cores gives 


/ = ja 2 Kp s d h 


(31) 


where re = irh/m is the circulation quantum, p s the superfluid density and the mutual friction parameter. Taking 
a from Fig. [2] and d|| from Ref. HI gives / = 1.6 x 10 -19 J s/m at the conditions stated above. 


IV. TWIST RIGIDITY 

The twisting of a double-core vortex is measured by a twist wave vector k = d z <f>. For a small twist the additional 
free energy is supposed to be quadratic in k, 


Etwist = \Kk 2 . (32) 

Here we aim to calculate the twist-rigidity coefficient I\. 

The reduced order parameter is expected to satisfy the boundary condition A (r, ip, z) —X \e tv for r —x oo, where I is 
a unit matrix. An axially symmetric vortex satisfies the following symmetry Q 

A(r, </) = e ie R(z, 0)A (r, 0 - 0)R(z, -0), (33) 

where R is a rotation matrix parametrized by an axis and angle of rotation. A twisted vortex apparently has the 
property 

A(r, (j), z) = e lkz R(z, kz)A(r , <j> — kz, 0)R(£, —kz). (34) 

Taking the z derivative and evaluating it at z = 0 gives 



( 35 ) 
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where 

^ 0-10 

10 0 

V° 0 0 

In the Ginzburg-Landau region the twist-rigidity coefficient K is given by 


~ ~ ~ f) A f) A 

B = iA + R'A — AR' + y^- — x—, R'= 
ox ay 


f 3 

3 

fd\J2 

J2\B m \ 2 + 2\B^ z \ 2 

m=i 

>=! 


(36) 


(37) 


In order to get a good approximation of the twist rigidity at general temperature, we have written the prefactor using 


Ag2 


h 2 p{ 1 - Y) 
40 mm* 


(38) 


This form of Aq 2 is valid since we assume Ff = Ff = 0. Physically, the twisting (1M1) causes axial spin flow in 
the asymptotic region and axial mass flow in the vortex center, but the Ff dependence of the former and the Ff 
dependence of the latter are not properly included in (1371) . 

The asymptotic form of the order parameter (neglecting dipole-dipole coupling) far from the vortex axis is of the 
form A = R(0) + 0(r~ 2 ), where 


9 = 


C\ cos 4> f sin </> 
1 + c 


r + cos 


C 2 sin 6 ( cos 6 

- 1 

1 + c 


sin (j)(j) 


(39) 


In general c = Aqi/2Ag 2 but in the present approximation Ff = Ff = 0 we have c = 1 [38|. Based on this, we can 
calculate the asymptotic contribution to the twist rigidity. Since (IMl) is expressed in cylindrical coordinates, the twist 
energy can be calculated more directly than above, by replacing 4>by (f>— kz. The additional energy can be evaluated 
from the gradient energy [3f| 


Fgz = 2 A G 2 (1 + c) J d 2 r^-^-= 2X G2 {1 +c)k 2 J d 2 r- 


= 2 Aq 2(1 + c)k 2 [ d(f) [ dr- 


dO k d9 k _ „ , c)k 2 [ ,00k 80 k 

+c)k J dr ww 

+ 4 cos 2 $ sin 2 (j) (Ci — C 2 ) 2 


dz dz 
(cos 2 cj) — sin 2 cj)) 2 


(1 + c) 


_ . 2 2 + 2c + c 2 . . 2 f j 1 

27rAG2fc --—-- Cl - C 2 / dr-. 

1 + c Jr 


Comparison to (13211 allows to identify 


Kt ail = 47TAG 2 2+ ^ + C \ ci - C 2 ) 2 [ dr 
1 + C J 


(40) 


(41) 


The dipole-dipole energy suppresses the contribution at distances beyond the dipole length 

We estimate numerical values under conditions explained above [below Eq. (1301) ]. Substituting the numerical order 
parameter into (1371) gives K = 3.2 x 10 -27 Jm in the region r < 60£o- Adding the asymptotic part from (l4lT) with 
Ci = 3.7-Ro and C 2 = 0.H? o we get K = 4.0 x 10 -27 Jm. 

Another estimation of the twist rigidity is to calculate the twisting energy of a pair of half-quantum vortices. 
Estimating energy by the length increase of the twisted pair gives 


IX p s n 2 a 2 u a 

-‘'■mass flow 397 f ^ f 


(42) 


where a is the distance between the half cores and r c their radius. Estimating a = 5.4I?o gives A" mass a ow = 1.9 x 10 -2 ' 
Jm. Estimating this without Fermi-liquid interaction, i.e. as we did in (1371) . we would get F mass fl ow o = 0-6 x 10~ 2 ' 
Jm. Supposing this difference is just missing in the formula based on (1371) . the corrected value for total twist coefficient 
is K = 5.3 x 10" 27 Jm. 
















